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Abstract 

We consider two simple models for the formation of polymers where at the initial time, 
each monomer has a certain number of potential links (called arms in the text) that are 
consumed when aggregations occur. Loosely speaking, this imposes restrictions on the 
number of aggregations. The dynamics of concentrations are governed by modifications of 
Smoluchowski's coagulation equations. Applying classical techniques based on generating 
functions, resolution of quasi-linear PDE's, and Lagrange inversion formula, we obtain 
explicit solutions to these non-linear systems of ODE's. We also discuss the asymptotic 
behavior of the solutions and point at some connexions with certain known solutions to 
Smoluchowski's coagulation equations with additive or multiplicative kernels. 
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1 Introduction 



The coagulation equations of Smoluchowski [T7] describe the evolution of concentrations of 
particles in a medium in which particles coalesce pairwise as time passes. The surveys by 
Aldous [T\ and Laurengot and Mischier [13] provide stimulating introductions to the subject. 
We also refer to Drake [6] and Dubovski [7] for detailed accounts, and e.g. to [9j [101 [15] and 
references therein for some more recent works in this area. 

In the discrete version of this model, particles are characterized by an integer m > 1 which 
should be though of as a size or a mass, in the sense that the result of the coagulation of a pair 
of particles, say {m, m'}, is m + m' . The dynamics are determined in terms of some symmetric 
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kernel k : N* x N* ^ K + such that, loosely speaking, every pair of particles {m,m'} coalesces 
at rate K(m,m'). This means that if c t (m) denotes the concentration of particles with size m 
at time t, then the evolution of concentrations is governed by the infinite system 



d 



m—1 



c t (m) = — c t (m')c t (m — vtl)k{vti \m — m') — c t (m) c t (m')K>(m,m') 



dt y ' 2 

m'=l m'=l 

where the first term in the right-hand side accounts for the creation of particles of size m by the 
coagulation of a pair of particles {m', m — m'}, the factor | stems from an obvious symmetry, 
and the second term accounts for the disappearance of particles with size m as the result of a 
coagulation with other particles. 

In general (OQ) cannot be solved explicitly; however there are a few important exceptions. The 
fundamental examples of solvable Smoluchowski's equations occur when the coagulation kernel 
k is constant, additive, or multiplicative. For instance, for a monodisperse initial condition 
(this means that at time t = 0, all particles have unit size), it is known since the original work 
of Smoluchowski [T7j that 

£\ -2 / j. \ m-1 



ct(m) = ( 1 + - 1 ( — j , < t < oo (2) 

solves ([1]) for K>(m,m') = 1. Further, if 

B{\, m) = (Am) m - 1 e- Am /m! , m G N* 
denotes the Borel probability function with parameter A G [0,1], then, according to Golovin 

c^m) =e-*B(l-e _ *,m), < t < oo (3) 
is a solution to for n(m, m') = m + m! . Finally, it has been shown by McLeod [H] that 

c t (m) = m~ l B{t, m) , < t < 1 (4) 

solves (CQ) for /t(m, m!) = mm'. More generally, explicit though complicated expressions can 
be obtained for solutions to flTJ when the coagulation kernel k is a linear combination of these 
three basic kernels; see [2UI HS1 EI] • 

The key reason why Smoluchowski's equations can be solved explicitly in these cases is that 
considering generating functions transforms ([T|) into a solvable PDE; see e.g. Deaconu and 
Tanre [5]. The purpose of the present work is to point out that the same techniques can be 
applied successfully to investigate a related but more complex model for coagulation. Roughly 
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speaking, the dynamics described by Smoluchowski's equations depend only on an additive 
quantity, namely the size of particles. Norris [16J developed the much more sophisticated 
setting of cluster coagulation models in which the rates of transition may depend on further 
parameters, such as typically the shapes of the clusters of particles. In this direction, we 
will deal here with a toy model for the formation of polymers in which coagulation rates are 
functions of a non-additive quantity that depends on the history of polymers. More precisely, 
the medium consists initially in monomers (elementary particles with unit size), such that each 
monomer has a certain number of potential links which we call arms. Specifically, we suppose 
that each particle is determined by a pair of integers (a, to), where a > is the number of arms 
of the particle and m > 1 represents its size. The arms serve to perform aggregations, i.e. to 
connect pairs of particles, and are consumed each time an aggregation occurs. We will consider 
two different models. 

The first model will be studied in Section 2. It is oriented, in the sense that a coagulation 
occurs when a particle uses one of its arms to grab another particle. Only one arm is consumed 
for each coagulation event; specifically, when a particle (a,m) with a > 1 grabs a particle 
(a', to'), these two particles merge into a single one (a + a' — 1, to + to'). We assume further that 
each arm grabs other particles at uniform rate; the precise formulation will be given in Section 
2.1. In Section 2.2 we use generating functions to connect these dynamics to a quasi-linear 
PDE that is then solved by the method of characteristics. Section 2.3 presents a version of 
Lagrange inversion formula which will be needed to invert the generating functions. The main 
results on the oriented model are given in Section 2.4, and finally some illustrative examples 
are discussed in Section 2.5. In particular, we point at some resemblances between on the one 
hand, the oriented model when the number of arms of monomers is distributed according to 
the standard Poisson law, and on the other hand, Smoluchowski's coagulation equation for the 
additive kernel. 

The second model will be studied in Section 3. It is symmetric, in the sense that each coag- 
ulation consumes one arm for each of the two particles involved into a coagulation event : 
when a pair of particles {(a, to), (a', m')} with aa' > 1 coalesces, the resulting particle is 
(a + a' — 2, to + to'). Furthermore, we assume that each pair of arms is activated uniformly. 
Section 3.1 introduces the setting rigorously, and points at a critical time in the evolution 
of the system which resembles the gelation time (i.e. the instant when particles with infinite 
size appear in the medium) in Smoluchowski's coagulation equations. We also use generating 
functions to make the connexion with a quasi-linear PDE closely related to the one that arises 
for the oriented model. The main results and some examples are given in Section 3.2. For 
instance, we determine the critical time in terms of the initial data (interestingly, r ra may 
be finite or infinite) and note some similarities between the symmetric model started from 
monomers with arms distributed according to the standard Poisson law and McLeod's solution 
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to Smoluchowski's coagulation equation for the multiplicative kernel at gelation time. 

We conclude this Introduction by explaning the term limited in the title. Roughly speaking, 
the effect of introducing arms in a coagulation model is that it imposes some restrictions to 
the number of aggregations. Indeed, in the oriented model, we shall show that if at the initial 
time the mean number of arms is less than the total concentration of particles, then the system 
converges as time tends to infinity to some limiting concentrations (and more precisely the 
system does not create particles with infinite size); see Corollary [TJ A similar phenomenon is 
observed for the symmetric model; see Corollary [2] and the remark thereafter. Perhaps this is 
more surprising for the symmetric model, as the latter bears some similarities with the classical 
Smoluchowski's coagulation for the multiplicative kernel, and it is well-known that gelation 
then always occurs in that case. 



2 The oriented model 
2.1 Setting 

We now describe precisely the dynamics for oriented coagulation with arms. For every t > 0, 
a G N := {0, 1,2,.. .} and m E W := {1, 2, . . .}, let c t (a, m) denote the density of particles 
(a, m) at time t. Loosely speaking, we suppose that at any time, each arm may grab each 
particle at the same unit rate. Rigorously, this means in terms of concentrations that the 
transition 

{(a, m), (a', m')} — ► (a + a' — 1, m + m) (5) 
occurs at time t with intensity 

ac t (a, m)c t (a', m) + a'c t (a, m')c t (a, m) = (a + a')c t (a, m)c t (a', m'). 

The evolution of the concentration functions is thus governed by an infinite non-linear system 
of ODE's 

^ j a+l m— 1 

—c t (a,m) = - / / (a + l)ct(a',m')ct(a — a' + l,m — m') 
at ' 2 ^-^ 

a'=0m'=l 

oo oo 

—c t (a, m) (a + a')c t (a', m!) . (6) 

a'=0 m'=l 
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We point out that 

a+1 m—l a+1 m—l 

2 



j a+1 m—l a+1 m—l 

- (a + 1)q(o/, m')c t (a — a' + 1, m — m') = a'c t (a', m')c t {a — a' + 1, m — m!) , 

a'=0m'=l a'=0m'=l 

and since in the sum in the right hand side, the term corresponding to a' = is null, we can 
reformulate as 

^ a+1 m—l oo oo 

— c f (a, m) = dct(a, m')ct{a — a' + 1, m — m!) — c t (a, m) (a + a')c t (a', m') . 

a'=l m'=l a'=0 m'=l 

The latter expression may be closer to the asymmetric description of the dynamics of the 
system. 

We also introduce the notation 

oo oo 

C t := ^2 c <( a ' m ) 

a=0 m=l 

for the total concentration of particles at time t, and 

oo oo 

A := ^ ^ ac t (a,m) 

a=0 m=l 

for the mean number of arms, and stress that (JH]) only makes sense as long as A t and C t are 
finite. But it is easy to check (see Lemma [T] below) that it suffices to assume that at the initial 
time Aq and Co are finite. 

It will be convenient to re-express the system ([6]) in a slightly different form by introducing 
the notation 

oo oo 

(°u f) ■= ^( a ' m ) c *( a ' m ) ' 

a=0 m=l 

where / : M x N* — > R+ stands for a generic nonnegative function. 

Lemma 1 For any solution to ([6]), the function t — > A t +Ct is non-increasing; as a consequence 
At + Ct < A + Cq /or allt>0. Further we have 



oo oo 



a,a'=0 m,m'=l 



whenever f : N x N* — > R zs bounded. 

Proof: First observe that (jH]) is precisely (j7j) when / is a Dirac function. By linearity, (j7j) 
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thus holds when / has finite support. 
Next fix n G N* and define 

if 1 + a < n, 
if n < 1 + a < 2n, 
otherwise. 

Then consider the function with finite support f n (a,m) := l{ TO < n j.<^ n (a). It is readily checked 
that 

f n (a + a' - 1, m + m') - f n (a, m) - f n (a', m!) < 

for all a, a' G M and m,m' G N*. We deduce from (ED) that t — > (q, f n ) is non-increasing. 
Further, f n (a,m) increases to 1 + a as n j oo, so by monotone convergence we have that 
A t + C t = lim n f tX) (c t , f n ), and we conclude that t — >• A t + C t is non- increasing. 

As we already know that ((Tj) holds when / has finite support, a standard approximation 
procedure yields the extension when / is only bounded by dominated convergence. □ 

We now start the analysis by pointing at an elementary property which should be intuitively 
obvious, as at each coagulation event, the number of particles and the total number of arms 
both decrease by one unit. The rigorous proof requires however some care. 

Lemma 2 For any solution to (jSJ), the difference D := Ct — At remains constant as time 
passes. More precisely : 

(i) If C = A , then 

C t = A t = C /{l + tC ). 

(ii) IfC -A = D^0, then 

r p Dt 

C - A ' + D = D C^-l ) + D - 

Proof: We get from ([7j) for / = 1 that 4iC t = —A t C t , and in particular the total concentration 
decreases with t. Similarly, fix a > and apply (J7j) for f^(a,m) = a A a to see that A^ : = 
(c t , /^) fulfills ^A^ < —A t C t . As At = lim^oo j A^ a \ this entails that A t also decreases 
with t. 

Then observe that 

| f a) (a + a' - 1, m + m') - f {a) {a, m) - f a) (a', m) \ < (a A a') + 1 , 
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and that therefore 

\f (a) (a + a'-l,m + rri) - f {a) (a, m) - f {a) (a', m) \ (a + a') < 2aa + a + a' . (8) 

Since 

lim {f {a) {a + a' - l,m + m') - f (a) (a, m) - f {a) {a ,m')) = -1 , 

a— >oo 

it is easy to get from ([7j) by dominated convergence that ^A t = —A t C t . We conclude that the 
difference C t — A t is constant and the explicit expressions in the statement follow readily. □ 

2.2 A quasi-linear PDE 

The system of evolution equations (j7j) resembles Smoluchowski's coagulation equation for the 
additive kernel. In the latter case, it is well-known that considering exponential functions yields 
a quasi-linear PDE related to the transport equation which can be solved explicitly, see e.g. 
[5]. This invites us to introduce the generating functions 

g t (x,y) := (c t , f x , y ) with f XjV : (a,m) -> x a y m and x,y G [0, 1] . 

For the sake of simplicity, we will focus on the situation where the total concentration of 
particles equals 1 at the initial time, i.e. Co = 1. The general case can easily be reduced to 
that one by a linear time-substitution. 

Lemma 3 (i) Consider initial concentrations (co{a,m) : a G N and m G N*) such that Co = 
Aq = 1. The system ([6]) is then equivalent to the quasi-linear PDE 

±!H(X,V) = (*(*,*) - ^9t(x,y) - j^g^y) 

for the generating functions, where t > andx, y G [0, 1]. The latter possesses a unique solution 
which is given in terms of its value at boundary t = by 

9t(x, y) = (1 + ty 1 g (h t (x ) y),y) = t~ [ h t (x, y) - , 

where h t (-,y) : [0, 1] — > [0, 1] is the (unique) inverse of the function x — > (1 + t)x — tgo(x,y). 

(ii) Consider initial concentrations (co(a,m) : a G N and m G N*) such that Cq = 1 and 
Co — A = D for some D G (— oo, 1)\{0}. The system (jBJ) zs then equivalent to the quasi-linear 
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PDE 

d , . ( . . De Dt \ d . . / De Dt \ , 

— qAx, y) = Qt(x, V) — x—pr- T^Qtix, y) + \ L> — Qt(x, y) 

dt K J V e - 1 + D J dx V e Dt -1 + D J y y ,y/ 

for the generating functions, where t > and x, y G [0, 1] . T/ie latter possesses a unique solution 
which is given in terms of its value at boundary t = by 

De Dt 

9t(x,y) = -~ — ——g (h t (x,y),y) 
e ut — 1 + D 

De Dt , , , D 2 e Dt 



e D '-l rv (e D '-l + J D)(e D *-l 

where ht(-,y) : [0, 1] — > [0, 1] zs i/ie (unique) inverse of the function 

x ^ D -i _ ! + D ) 3. _ _ x ) ^ y) ) 

Proof: (i) From ((7j) for / = / x>y , we get 

d 



(it 9t( x ,y) = {a-/3-j)/2 



where a is given by 



oo oo oo oo 



CLX a - 1 y m c t {a,m)x a 'y m 'c t {a f ,m')+ ^'"V' ' c t {a! ,m)x a y m c t {a,m) 



a,a—0m,m=l u,o=0m,tn=l 



d 

2gt(x,y)—g t (x,y) 



and /? = 7 by 



oo oo oo oo 

X 

a,a'=0 m,m'=l a,a'=0m,m'=l 



^ 2J ax a 1 y m c t (a,m)c t (a',m') + /J /J a'x a y rn c t (a,m)c t (a l \w!) 



d 

= xCt—gt{x,y)+gt(x,y)A t . 

Recall from Lemma [2] that C t = Af = 1/(1 + 1) to get the PDE in the statement. 

Note that this PDE does not involve partial derivatives with respect to the variable y, and 
is quasi-linear. It can be solved using the methods of characteristics. Specifically we consider 
functions t = t(r, s), x = x(r, s) and u = u(r, s) for < r < 1 and s > 0, with the boundary 
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conditions t(r, 0) = 0, x(r, 0) = r, u(r, 0) = go(r, y) and such that 

dt dx x du u 

1 — = u 



ds ' ds 1 + 1 ' ds 1 + 1 

We deduce 

t(r,s) = s , u(r, s) = g (r,y)/(l + s) , x(r, s) = r(l + s) - s^ (r, y) ; 

and then the comparison with the quasi-linear PDE gives gt(x, y) — u(r,s). 

The requirements C = g (l,l) = I and A = -^g (l, 1) = 1 ensure that for any s > and 
y G [0, 1], the function r — > r(l + s) — sg (r,y) has derivative 1 + s(l — -^go(r^y)) > 1, and 
since — sgo(0, < < 1 < (1 + s) — sgo(l, y), there is a unique solution h s (-, y) : [0, 1] — > [0, 1] 
to the equation h s (x, y) = h s (x, y)(l + s) — sg (h s (x, y),y). This provides the solution given in 
the statement. 



Conversely, if (gt ■ t> 0) is a family of generating functions solving this PDE with a boundary 

dx ' 



condition go(x, y) = (cq, f XiV ) and such that go(l, 1) = Cq = 1 and -§-go(l, 1) = A = 1, then we 



know from above that 

x 



9t(x, y) = (1 + 1) 1 g {h t (x,y),y) =t 1 h t (x,y) 



t 2 + t 



with ht(-, y) : [0, 1] — ► [0, 1] the (unique) inverse of x — > tgo(x, y). Note that 1) = 1 

for all t, which gives g t (l, 1) = 1/(1 +t). Further, taking a (left) derivative with respect to the 
variable x at x = 1 in the identity 



x = {l + t)ht(x,l)-tg Q (h t {x,l),l) 



entails J^h t (l, 1) = 1 and then 

^9t(l, 1) = r 1 ^ 1 ' !) " V(t 2 + *) = V(l + • 
Thus the family (g t : t > 0) also solves the equation 

d d d 

fa9t(%, y) = (gt(x, y) - mi 1 , 1)) g^9t(x, y) - g t (x, y)—g t (l, 1) . 

Since ^(1,1) = (c t , 1) = C t and ^^(1,1) = (c t , a) = A t , we recover (jSJ) by inverting the 
generating functions. 

(ii) The derivation of the PDE is obtained as in the proof of (i), applying the second part 
of Lemma [2] in place of the first. It is solved again by the methods of characteristics; using the 
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same notation as in (i), we are led to consider the system of ODE's 

dt dx _ De Dt du _ f e Dt \ 

d^" 1 ' te~ X e Dt_ 1 + D ~ u > ^ ~ V ~ e D * - 1 + D ) U 

with the boundary conditions t(r, 0) = 0, x(r, 0) = r, u(r, 0) = go(r, y). We obtain 

De Ds 

t(r,s) = 8 , u(r, s) =9o(r,y) eDs _ 1 + D 

and 

x(r, s) = D- 1 (r (e Ds - 1 + D) - g Q (r, y) (e Ds - l)) . 
The proof can then be completed just as in the case (i). □ 

2.3 A version of Lagrange inversion formula 

The final step of our analysis consists in checking that the function (x,y) — > h t (x,y) that 
appears in Lemma [3] is the generating function of some finite measure on N x N* that can be 
inverted explicitly, at least under some natural hypotheses. It is an easy application of Lagrange 
inversion formula (see, e.g. Section 5.1 in [23]). 

Lemma 4 Let fi = (//(a)) a6 N be a finite measure on N, \i ^ 0. We write 

oo 

g(x) = ^x a /i(a) , x G [0, 1] 

a=0 

for its generating function. For every p, q > and for every x, y > sufficiently small, the 
equation 

h(x,y) = yg(px + qh(x,y)) 
has a unique solution which is analytic in x and y and given by 

OO OO 1 / I _ i \ 

h(x, y) = J2J2 x<1 y m - g m ~ V/i* m (m + a - 1) , 

where fi* m stands for the m-th convolution power of [i. 

Proof: We fix p, q > and x sufficiently small so that px < 1. Then we define g(y) : = 
g(px + qy) for y > with px + qy < 1, which is the generating function of the sigma-finite 
measure 

OO / \ 

n 



n=k 
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Observe also that g(0) = g(px) > 0. According to the Lagrange inversion formula, the equation 
h(x,y) = yg(h(x,y)) has a unique solution for y > sufficiently small, which is analytic in y 
and can be expressed as 

oo 

h{x,y) = Y,y n n- 1 f n {n-l). 

n=l 

Then observe that the generating function of the measure fi* n is g n (y) = g n {px + qy), so 



iT(k) = E ( ; ) w~V/i* n (j) 

j=k 



We deduce that 



h(x, y ) = e^- 1 Yl ( n J _ 1 ) (^r n+ v~y n (j) 



j 

\ n 

n=l j=n—l 

and the change of variables j = a + n — 1 completes the proof. □ 



2.4 Explicit solutions 

We are now able to solve ([6]) explicitly in the situation when the initial concentrations co(a,m) 
are carried by particles having unit size; recall also the assumption that the initial total concen- 
tration is 1. In this direction, it may be interesting to point out that even though the number 
of arms of a particle is not an additive quantity, the simple change of variables a = a — 1 yields 
a new parametrization of particles that is additive on coagulation, in the sense that ([5]) then 
reads 

{(a, m), (a', m)} — > (a + a', m + m') . (9) 

This observation makes the solvability of ([6]) by techniques based on generating functions easier 
to understand. Note however that the parameter a may take the value —1 and that the rate 
of the coagulation is now a + a' = a + oi + 2 and hence not additive in the alternative 
parameter a. 

More precisely, we consider ji = (/i(a)) ae N a probability measure on N, and denote for every 
m G N* by /i* m the m-th convolution power of \i. 

Theorem 1 The system ([6]) has a unique solution (c t (a,m) :oeN,m6f and t > 0) started 
from 

c (a,m) = l{ m=1 }/x(a) , a G N and m G N* 
which is given as follows : 
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c t (a,m) = m-H m - 1 (l + t)- (a+m) ( ° + ™ 1 ] // m (a + m - 1] 



(ii) If J2aen a M a ) = I — D for some D E (—00, 1)\{0}, th 



en 



,,,.„„ = (e» - 1)- (e<* - 1 + D)™ ( « + <" " 1 ] + m - 1) 



Proof: We shall focus on (i), the argument for (ii) being similar (note also that taking the 
limit as D — > in part (ii) yields the solution in the case (i), as it should be expected). 

According to Lemma [3] and the present hypotheses, we have to consider the initial generating 
function 

oo 

9o(x,y) = (co,fx, v ) = ^x a yn{a) = yg(x) 

a=0 

where g is as in Lemma HI and then the solution h t (x, y) to 

(1 + t)ht(x, y) - tyg(h t (x, y))=x. 
Lemma H] invites us to introduce 

h(x,y) := t _1 ((l + t)h t (x,y) - x) , 

so that 

h(x, y) = yg (^~[ + Y+t V >) ' 

According to Lemma 0] (with p = 1/(1 + 1) and q = t/(l +£)), the solution to this equation can 
be expressed in the form 

^») = EE^ [ a (rh) (i+rv m (m+a-i), 

whenever x, y are sufficiently small. Recall that 

h t (x,y) = (l + t)- 1 (th(x,y)+x) : 
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so we get 



00 00 1 / 1 1 \ 

^(z, y ) = + V V 2r a 2/ m - m a ~ t m (l + t)-( a+m V* m (m + a - 1) 

1 + * „ , m \ a / 

a=0 m=l \ / 

We then known that the unique solution to the quasi-linear PDE in Lemma [3](i) is 

-l? ' x 



gt(x,y) = t h t (x,y)- 



= VViT- | ^(l + ^-^+^u^fm + a-l). 

^ ^ m \ a / 

a=0 m=l \ / 

Hence g t coincides with the generating function of the concentrations that appear in the state- 
ment, and the proof can be completed by an appeal to Lemma E](i). □ 

We now conclude this section describing the limiting behavior of concentrations as time 
tends to infinity, starting with the case D > (i.e. at the initial time, the mean number of 
arms is less than the total concentration). 

Corollary 1 Let // = (/i(a)) a eN be a probability measure on N with mean 

y a/i(a) = 1 - D < 1 . 

Then the solution (cj(a, m) : a G N, m G N* and t > 0) to the system started from 

Co (a, m) = l{ m= i}/z(a) , a G N and m G N* 
/ias a limit as t —>■ 00 raf^NxN*) gwen fry 

Coo(a,m) = l {a=0} — /i* m (m - 1) . 
m 

Proof: The pointwise convergence (i.e. with a and m fixed) should be plain from the expres- 
sion given in Theorem [Tfii). Recall from Lemma [2] that the total concentration C t tends to D 
as t — > 00. On the other hand, as /i is a probability measure on N with mean less than 1, it is 
easily seen (see e.g. the remark after this proof) that m _1 /i* m (m — 1) also defines a probability 
measure on N*. So 

00 00 

a=0 m=l 

and we can complete the proof invoking Scheffe's lemma (see for example [1]). □ 
It is well-known and easy to check that if jj, is a probability measure on N with mean 
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1 — D < 1, then \x := {m l n* m {m — l)) mG N* is a probability measure on N* with mean 

oo 

J>*™(m - 1) = LT 1 . 

m=l 

Indeed, we know from the Lagrange inversion formula (see for instance Theorem 5.1.1 in |23j) 
that the generating functions g of /i and h of pL are related by the equation h(x) = xg(h(x)). 
The condition on the mean of fi reads g'(l) = 1 — D < 1; this readily entails that h(l) = 1, so 
that p, has total mass 1. Then taking the (left) derivative at x — 1, we get ft/(l) = g(h(l)) + 
ft/(l)</(/i(l)), which yields /i'(l) = 1/D. Thus Corollary [1] implies that 

oo 

^ mcoo(0, m) = 1 . 

m=l 

This property can be interpreted as follows. Suppose that at the initial time, we tag a monomer 
uniformly at random. Then the distribution of the polymer at time t that contains this tagged 
particle has the distribution 

oo oo 

} u ^ mc t (a, m)6(a, m ) , 

a=0 m=l 

and as time passes, this family of measures remains tight. Physically, this means that if at the 
initial time the mean number of arms is less than the total concentration of monomers, then 
the oriented model does not produce particles with infinite size as time tends to infinity. 

Remark. The recent paper [3] deals with a system of randomly interacting particles which is 
closely related to the present deterministic model, and sheds a probabilistic light on Corollary 
[TJ More precisely, in [3], time is discrete and at the initial time there are n particles with 
arms such that the sequence £i, . . . , £ n of the number of arms of particles is i.i.d. with a fixed 
distribution /x on N. Arms are enumerated uniformly at random, which specifies the order of 
activation. When an arm is activated, it grabs uniformly at random one of the particles in the 
system which had not been grabbed previously and which does not belong to its own cluster 
either. The polymerization procedure terminates when all arms have been activated, and the 
terminal configuration is given by a forest of trees. Roughly speaking, the main result in [3] is 
that if /i is subcritical (i.e. its first moment is 1 — D < 1), then as n -» oo, the distribution of 
a tree picked uniformly at random in the terminal configuration converges to that of a Galton- 
Watson tree with reproduction law //. According to Dwass [8], the probability that the size of 
a Galton- Watson tree with reproduction law // is m equals m _1 /i* m (m — 1). Further, by the 
law of large numbers, the number of trees in the terminal forest is approximately Dn, so that 
the density of trees with size m at the terminal time is Dm~ l /i* m (m — 1). This corroborates 
Corollary [TJ proving another example of the deep connexions between coagulation models and 
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branching processes (see, e.g. [5| [19]). 



2.5 Examples 

Let us now discuss some explicit examples. 

Consider first the degenerate case when \i = S± is the Dirac mass at 1, so at the initial time, 
there is a unit concentration of particles with unit size having exactly one arm, and all the 
other concentrations are 0. It is clear that coagulations then always produce polymers with 
exactly one arm, and therefore in this specific situation, the notion of arms plays no role in the 
evolution of the system. In particular it is not surprising that a known solution should emerge, 
however it may be interesting to discuss this case as a verification of our general formulas. 
Specifically we have /x* m = 5 m so we find 



Q (i,m) = r- 1 (i + t) 



and c t (a,m) = for a / 1. As a check, observe that summing these quantities for m G N* 
is in agreement with Lemma [3J It is also interesting to compare with Smoluchowski's solution 
f$2§. Specifically we see that Smoluchowski's coagulation equation for the kernel K(m,m') = 2 
and monodisperse initial condition can be viewed as the present oriented coagulation with arms 
when at the initial time all particles have unit size and exactly one arm. This observation can 
be established directly by an elementary analysis of the transition rates in both models. 

More generally, consider the case when jj, is the binomial law with parameter (n, 1/n), where 
n > 2 is some integer, i.e. 

^(o) = rM n' n (n - l) n ~ a for a = 0, . . . , n . 

Then fi* m is the binomial distribution with parameter (mn, 1/n), and we get 

ct(a, m) = ^— ^ — — |-rt m_1 (l + t)- {a+m) n- mn (n - 

[ran + 1 — a — m)\m\a\ 

when a + m — 1 < mn, and c t (a, m) = otherwise. 

Then we let n tend to oo and thus consider the case when fi is the standard Poisson law, i.e. 
n(a) = l/(a!e) for a G N. Then /i* m is the Poisson distribution with parameter m and we get 

c t {a,m) = e- m m a+m - l t m - l (l + t)- (a+m) , m G N. 

a\m\ 

It is interesting to point out that summing this quantity over a G N yields the total concentra- 
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tion of particles with size m, 

C t {m) := 5>(a,m) = -J- \ ^^ e^) = ^-5(^/(1 + t),m) . 

The comparison with Golovin's solution ([3]) suggests that, loosely speaking, Smoluchowski's 
coagulation equation for the additive kernel and monodisperse initial condition coincides after 
the logarithmic time change t — > log(l + t), with the present oriented coagulation with arms 
when at the initial time, particles have unit size and the number of arms is distributed according 
to the Poisson law. 

We leave to the interested reader the task of developing similar calculations for binomial or 
Poisson laws with mean ^ 1 to illustrate Theorem [T^ii). Here is a final example in this vein. 
We assume that p, is the Negative Binomial distribution with parameters r > and p G (0, 1), 
viz. 

fl ( a ) = 1 %±^p r (i-py, a en. 

all (r) 

Recall that 

oo 

= r ( l - P)/P :=1-D, 

a=0 

and that p* m is Negative Binomial distribution with parameters mr and p. Assuming that 
D 7^ 0, we then find that the concentration c t (a, m) is given by 

e Dt m -i D a + i uot _ jxm-i urn _ j + ^ -(»+») ?{mr + a + m - 1) _ _ x 

v ' v ' a\[m — ly.T(rm) 



3 The symmetric model 

3.1 Setting and relation to a quasi-linear PDE 

We next turn our attention to the symmetric model of coagulation with arms, keeping the 
notation for the oriented one. This means that now each aggregation event consumes two 
arms, one for each particle involved, and that any pair of arms is activated at the same unit 
rate. So the transition 

{(a, m), (a', m')} — ► (a + a' — 2, m + m ) (10) 
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occurs at time t with intensity 

aa'c t (a, m)c t (a', m') , 

and the evolution of the concentration functions is specified by the infinite non-linear system 
of ODE's 



— Ct(a, m) = — a' {a — a' + 2)c t (a f , m')c t {a — a' + 2, m — m' 



a'=l m'=l 



oo oo 



c t (a,m)2_^ aa ' c t( a> , m ')- (11} 



a'=l m'=l 



Again (TTTj) only makes sense as long A t < oo, but it is readily seen that t ^ A t decreases as t 
grows, so it suffices to require that A < oo. 

Before starting the analysis, let us point at the special role of particles with no arms. Indeed, 
particles with no arms are inactive in the symmetric model, in the sense that they cannot 
coagulate with other particles (this was not the case in the oriented model as a particle with 
no arm could still be grabbed by some other particle). Analytically, this is seen from the fact 
that the sub-system ffTTj) for (a, m) G N* x N* is autonomous. We also stress that particles with 
no arms are produced by the coagulation of two particles both with a single arm, and more 
precisely, specializing ( TTTj) yields the simple identity 

d 1 m 1 

— q(0, m) = - 22 °t{^ m ') c t( 1 ^ m-m!) . (12) 

m'=l 

Just as for the oriented model, it is convenient to re-express the system ( TTTT) as 

^ ^ oo oo 

^( c ^/) = 2 (f( a + a '- 2 , m + m ')-f( a , m )-f( a '^ m l) aa '^( a , m )^( a ', m '), ( 13 ) 

a,a'=l m,m'=l 

where / : N x N* -> R + is a generic nonnegative and bounded function. 

The equations (TTTj) resembles Smoluchowski's coagulation equation ([1]) for the multiplicative 
kernel. In the latter case, it is well-known that a phenomenon of gelation occurs, in the sense 
that the total mass is not a preserved quantity for all times as one might expect naively. 
Informally, this is due to the formation of particles of infinite size in finite time; see e.g. 
[H[T2l[22]. A similiar phenomenon may (or may not) happen in the present case, and we shall 
study ( TTTj) and ( TT3j) before that critical time. Specifically, we introduce for any r > 

r r := inf {t > : (q, a 2 ) > r} , 
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where by a slight abuse in notation, we write a 2 for the function (a, to) — ► a 2 , and then 

:= sup{r r : r > 0} . 
Recall also that A t := (q, a) is the mean number of arms at time t. 

Lemma 5 For the symmetric model, we have 

A t = for all t < Too . 

l + tA J 

Proof: The argument is similar to that in Lemma [2j Using the same notation as there, we 
first specify (Tl3|) to the function and let a — > oo. The bound 

| (a + a' - 2, to + to') - (a, to) - (a', to') | < (o A a') + 2 , 

enables us to apply the theorem of dominated convergence provided that sup 0<s<t (c s , a 2 ) < oo 
and we get the equation ^A t = —A 2 . □ 

We stress that the formula in Lemma may fail when t is too large. Indeed, applying (1131) 
to f(a,m) = 2, we also get ^(2C t ) = —A 2 without requiring that t < Too. Thus if Lemma [5] 
was always valid, then the difference 2C t — A t would remain constant. But this is absurd when 
2C < A since then one would have 2C t = 2C + A t — A — » 2C — A < as t — > oo. 

For the sake of simplicity, we shall focus in the rest of this section on the case when A : = 
(c , a) = 1, which induces no significant loss of generality as the general case can be reduced 
to that one by a linear time-change (provided that of course A < oo). Recall the notation 
g t (x,y) for the generating function of c t and introduce 

d 

kt{x,y) := Q^9t{x,y) , x,ye[0,l] 
which should be viewed as the generating function : 

oo oo 

k t (x, y) = J2J2 xa y m ( a + 1 ) c *(° + m ) • 

a=0 m=l 

The following statement is a partial counter-part of Lemma [3] for the symmetric model. 

Lemma 6 Assume that A = 1. Then for any solution to the system (fTT!) we /iai>e the equation 

d 9 

— kt(x, y) = (k t (x, y) - xA t ) g^ k t(x, y) - A t k t (x, y) , 
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and for t < T^, the latter can be rewritten as the quasi-linear PDE 

^-k t (x,y) = \k t (x,y) — ] -^-k t (x,y) k t (x,y). 

dt v ,yj V 1 + t J dx K ,UJ 1 + t K ,yj 

Proof: Just as in the proof of Lemma [3], the equation f|T3|) specialized to f(a, m) = f x , y (a, m) = 
x a y m yields 

Then taking the partial derivative with respect to x, we obtain the first equation in the state- 
ment. The second follows from an application of Lemma [5] and the assumption Aq = 1. □ 

We observe that the second PDE in Lemma [6] is the same as in Lemma [3](i) for the oriented 
model with k t (x, y) replacing gt{x, y). In this direction, we recall that a similar relation between 
solutions of Smoluchowski's coagulation equations for the additive and the multiplicative kernels 
holds, see e.g. Theorem 3.9 in [5]. 



3.2 Explicit solutions and examples 

We are now able to solve (fTTI) up-to the critical time Too when at the initial time, all particles 
are monomers, i.e. each particle has unit size and its number of arms is arbitrary. Just as for 
the oriented model, it may be interesting to observe that the re-parameterization a = a — 2 is 
additive on coagulation, in the sense that (fTUl) then reads 

{(a,m),(a',m')} — > (a + a', m + m) (14) 

This also makes the solvability of ([6]) easier to understand; note however that the parameter a 
may take negative values and that the rate of the coagulation (fT4"l) is now aa' = (a + 2) (a' + 2) 
and hence not multiplicative in the alternative parameter a. 

We consider a measure fi = (/i(a)) a eN on N with unit mean and finite second moment (we 
stress that we do not require /x to be a probability measure), and introduce the probability 
measure v = (v(a)) aE N given by 

u(a) = (a + l)/i(a + l), a G N . 

We denote the first moment of v by 

oo oo 

M := av[a) = a(a — l)/i(a) , 

a=0 a=l 
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and then define 

f oo if M < 1 

T = < 

|1/(M-1) otherwise. 
Recall also that we write u* m for the m-th convolution product of za 

Theorem 2 T7ie system ffTTj) /ias a unique solution (c t (a,m) : a, m G N* and £ < T) started 
from 

co(a, m) = l{ m= i}//(a) , a, m G N* 
which is given for a, m > 1 fry 



Q (a,m) = v — , y " £ m - 1 (l+£)-( a+m ^ 1 )z/* m (a + m-2). 
a!m! 



(a + m 2)! ^ m _ 1 ^ ^ n _ l „ + ,„_ j , ;/ ,,„ f 
a!m! 

Further, T coincides with the critical time T^. 

Proof: The initial conditions invite us to introduce the generating function of /x, 

oo 

g(x):=^2x a fi(a), x G [0, 1] , 

a=0 

and to set for y G [0, 1] 

fltoOc, 2/) = yg{%) and fc (x, y) := q^9o(x, y) = y ^M )- 

a=0 

Consider the second PDE of Lemma El viz. 

^h(x, y) = (k t (x, y) - J~^*( x > V) ~ i~\~t^ Xj ^ ' ^ 

with boundary value ko(x, y) defined above. If we replace k t by g t , this is precisely the PDE that 
has been solved by the method of characteristics in Lemma E^i). There is however a difference 
that requires some attention : here we assume that -^k (l,l) = M < oo whereas we had the 
stronger hypothesis S,go(l, 1) = 1 in Lemma [3](i). Nonetheless, the condition £ < T implies 
1 + £ — £M > 0, and thus the derivative of the function x — > (1 + t)x — tko(x,y) is strictly 
positive for x G [0, 1]. As — tk (0, y)<0<l<l+£ — tk (l, y), this ensures the existence of a 
unique inverse function £ t (-, y) : [0, 1] — > [0, 1] to x — > (1 + t)x — tk (x, y). The argument in the 
proof of Lemma Efi) is thus still valid, and we conclude that for every t < T and x, y G [0, 1] 



fct(a;, y) = (l + t)- l k (£t(x, y),y) = r l i t {x, y) - 



t 2 + t 
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The comparison with Theorem []Ji) now yields that for t < T and x, y G [0, 1] 

OO OO / I _ 1 \ 

k t (x, y) = x a y m rrrH m - l {l + t)~ ia+m) i a + m „*m( a + m - 1) . 



a=0 m=l 

We thus see that if we define c t (a, m) for a, m G N* and < t < T as in the statement, and let 
kt be the generating function of ((a + + l,m) : a G N and m G N*), then fc t solves ( [151) 

for < t < T. 

We next observe that £ t (l, 1) = 1 since k (l, 1) = 1, and thus 

k t (l,l) = (c u a) = l/(l + t). 

Hence k t also solves 

d .9 

— k t (x, y) = (k t (x, y) - x(c t , a)) Q^ k t( x , v) ~ ( c t, a)k t (x, y) , 

which is the first PDE in Lemma [6j Inverting the generating functions, we conclude that 
(ct(a, m) : a,m G N* and < t < T) is a solution to ( TlTT) . 

We then check that the critical time of this solution coincides with T. In this direction, 
we first recall that -^k (l,l) = M, £ t (l,l) = 1 and, by definition, that 

{l+t)£t(x,l)=tko(£t{x,l),l) + x. 

We take the (left) derivative with respect to the variable x at x — 1 and obtain 

(i + t)— e t (i,i) = tM— ^(i,i) + i, 

ox ox 

so 

It follows that 

9x v ' ; 9x v ' ; v ; (l + *)(l + t(l-M)) 

remains bounded on compact intervals in [0, T[, and further explodes as t | T when M > 1 
(i.e. T < oo). As ^/c t (l, 1) = (Q,a 2 — a), we conclude that T = Too. Finally, the uniqueness 
of the solution to (1TT1) up to the critical time should be plain from a perusal of the preceding 
arguments. □ 

Theorem [2] does not provide an expression for concentrations of particles with no arms. 
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However, these can be recovered from the concentrations of particles with exactly one arm via 
the equation ()12p (note that the concentration of monomers with no arm does not evolve, so 
we may focus on particles with size m > 2 in the following statement). 

Corollary 2 Under the same assumptions as in Theorem^ and assuming also that z/(0) > 0, 
we have for every m > 2 and t < T 

c t (0,m) = — - -(1 + l/t) 1 - m u* m (m-2) form>2. 

m{m — 1) 

As a consequence, in the case when M < 1 (that is when T = oo), the concentration of 
particles (a, m) at time t has a limit c^a, m) as t — > oo which is given by 

Coo(a, m) = l|a=o> — 7 -u* m (m — 2) for a G N and m > 2 . 

m(m — 1) 

Proof: We know from Theorem [2] that for every t < T and m > 1 

q(1, m) = t m -\l + ty m m- l v* m (m - 1) . 

According to the classical Lagrange inversion formula (cf. Theorem 5.1.1 in [23J), m~ 1 z/* m (m — 
1) appears as the m-th coefficient in the analytic expansion of the entire function u which 
solves u(x) = x(p(u(x)), where is the generating function of the probability measure v (the 
assumption u(0) > ensures that (f)(0) > 0). Another application of the Lagrange inversion 
formula now shows that the m-th coefficient in the analytic expansion of the entire function 
u 2 (x) is 2m~ 1 v* m (m — 2). In terms of convolution product, this reads 

m-l 2 



q(1, -)* 2 (m) := V q(1, m'Wl, m - m!) = t m - 2 (l + t)- m — v* m (m - 2) . 
^— ' m 

m'=l 

Since Co(l, m) = for m > 2 by assumption and 

/ s m ~ 2 {l + s)- m ds = {m- + l/t) 1 -™ , 

Jo 



our first claim follows from f|T2l . The second follows immediately from the first for a = 0, and 
from Theorem [2] for a > 1. □ 

We stress that when M < 1 and v ^ 5%, then 

oo oo oo 1 oo 

E mC ~( a > m ) = E 7 -rv* m {m - 2) = E — rK«) = / / (^*) • (16) 

' ' (m — 1) a + 1 
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Just as for the oriented model, this identity can be interpreted as a tightness property as time 
passes for the distribution of the polymer at time t that contains a randomly tagged monomer. 
Physically, this means that the symmetric coagulation model does not produce particles with 
infinite size as time tends to oo. 

Indeed, the identity (ITB|) is easily checked as follows : 

OO oo oo 

V- u* m (m-2) = VK«)52t -^ m - l \m-2-a) 

t-^im—l) ^ ^ (m-l) v ' 

m=2 v > a=0 m=2 v ' 

oo oo 

= ^v(a)^n- l v* n {n-a-l) . 

a=0 n=l 



As v 7^ Si is a probability measure with mean M < 1, it follows easily from the Lagrange 
inversion formula (see the argument in the remark after the proof of Corollary [1]) that for any 
a > 0, 

oo 

^i±V»( n _ a _l) = l. 

n=l 

We stress that this is where the assumption v ^ Si is needed as otherwise the generating 
function of v would vanish at 0, impeding the application of Lagrange inversion formula. We 
thus have checked that 

oo 

V -v* m (m - 2) 

^ m-1 V ; 

m=2 v 1 

the other equalities in (Tl6]) are obvious. 

To summarize this observation, if we exclude the case when v — Si (or equivalently, fi = 
that will discussed further below, we have shown that there either gelation occurs at a finite 
time (by Theorem EJ this happens if and only if M > 1), or the mass distribution of particles 
remains tight as time tends to infinity. 

Remark. The limiting concentrations in Corollary [2] bear a striking resemblance with the dis- 
tribution of the total population generated by a (sub)-critical Galton- Watson branching process 
with reproduction law v and started from two ancestors. We refer the interested reader to the 
recent work [2] for a probabilistic interpretation that relies on the study of the random config- 
uration model. The latter can be used to create typical graphs with specified degree sequence, 
and is constructed via a stochastic algorithm which can be thought of as the probabilistic 
counterpart of the present symmetric model of coagulations with limited aggregations. 

We now conclude this work by illustrating Theorem [2] with some examples. The simplest is 
when fi = Si, that is, initially, there is a unit concentration of monomers with a single arms. 



00 

= 5> + l)-Ma), 

a=0 
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So v = 5q and one finds 

^(1,1) = ^ , c *(°> 2 ) = ^^' 0<t<T = oc, 

and all the other concentrations are zero. 

Next, suppose that at the initial time, there is a concentration 1/2 of monomers with two 
arms, i.e. /i = \8 2 - Then v = 8\ and we get 

ct(2, m) = -t m -\l + t)- {m+1) , meN*,0<t<T = oo, 

and all the other concentrations are zero. Note the similarity with the oriented model started 
from a unit density of monomers with a single arm. Of course, this property could also be 
observed by a direct argument. 

We then turn our attention to the simplest example with a finite critical time, namely the 
situation where at the initial time, there is a concentration 1/3 of monomers with three arms, 
i.e. /i = |$3. Then v = 82 and M = 2, so the critical time is T = 1. One gets 

c t (m + 2, m) = ( 2m ) ! f"-i(i + £)~(2m+i) m e N* and t < 1 , 
(m + 2)!m! 

and all the other concentrations are 0. It is easily checked that 

00 

lim ^^(m + 2) 2 ct{m + 2, m) = 00 , 

m=l 

in agreement with Lemma 

Finally, consider the case when \x is the standard Poisson law, then v is also the standard 
Poisson law and in particular T = 00. As v* m is the Poisson law with parameter m, we get 
from Theorem [2] that 

ct(a, m) = — r-^^ 1 ! 1 + ^-^-^e"" 1 , a, m e N*, t > , 
a!m! 

and from Corollary [2] that 

q(0, 1) = e" 1 and q(0, m) = e- m ^-^(l + lit) 1 "™ for m > 2 . 

m! 

It is interesting to point out that 

lim c t (0,m) = e" m — = m~ l B(l, m) , meN* 

t^oo ml 
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and to compare with McLeod's solution (j4j) to Smoluchowski's coagulation equation for the 
multiplicative kernel and monodisperse initial condition. We see that the terminal state for 
the symmetric model of coagulation started from monomers with arms distributed according 
to the standard Poisson law is the same as the state at gelation time in Smoluchowski's model 
for the multiplicative kernel and monodisperse initial condition. 
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